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Abstract 

We study the approximate analytical solutions of the Dirac equation for the generalized Woods- 
Saxon potential with the pseudo-centrifugal term. In the framework of the spin and pseudospin 
symmetry concept, the approximately analytical bound state energy eigenvalues and the corre- 
sponding upper- and lower-spinor components of the two Dirac particles are obtained, in closed 



qh' form, by means of the Nikiforov-Uvarov method which is based on solving the second-order linear 
differential equation by reducing it to a generalized equation of hypergeometric type. The special 

\ cases K = ±1 (/ = Z = 0, s-wave) and the non-relativistic limit can be reached easily and directly 
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I. INTRODUCTION 



It is well known that the exact solutions play an important role in quantum mechanics 
since they contain all the necessary information regarding the quantum model under study. 
The exact solutions of the Schrodinger equation are only possible for a hydrogen atom 
and for a harmonic oscillator in three dimensions [1-3]. However, when a particle is in a 
strong potential field, the relativistic effect must be considered, which gives the correction 
for non-relativistic quantum mechanics [4]. Taking the relativistic effects into account, a 
particle in a potential field should be described with the Klein-Gordon (KG) and Dirac 
equations. In recent years, there has been an increased interest in finding exact solutions 
to Schrodinger, KG, Dirac and Salpeter equations for various potential schemes [4-20]. 
The problems that can be exactly solved for the KG and/or Dirac equations are seldom 
except a few examples, such as hydrogen atom and electrons in a uniform magnetic field. 
Recently some authors solved such relativistic equations for some potentials. Ikhdair [6] 
obtained the bound-state solution of the L>-dimensional KG equation for the vector and 
scalar general Hulthen-type potentials with any arbitrary Z-state using the Nikiforov-Uvarov 
(NU) method. Moreover, Egrifes and Sever [7] investigated the bound state solutions of 
the ID Dirac equation with PT-symmetric real and complex forms of generalized Hulthen 
potential. Yi et al. [10] obtained the energy equations in the KG theory with equally mixed 
vector and scalar Rosen-Morse-type potentials. We have solved the spinless ID Salpeter 
equation analytically for its exact bound state spectra and wave functions with real and 
complex forms of the T'T-symmetric generalized Hulthen potential [11]. We obtained a 
quasi-exact analytic bound-state solutions within the framework of the position-dependent 
effective mass KG equation for scalar and vector Hulthen potentials in any arbitrary D- 
dimension and with any orbital quantum number / using the NU method combined with a 
new approximation scheme for the centrifugal potential term [20]. 

The Woods-Saxon (WS) potential and it's various modifications have played an essential 
role in microscopic physics in the determination of the energy level spacing, particle number 
dependence of energy quantities and universal properties electron distributions in atoms, 
nuclei and atomic clusters since it can be used to describe the interaction of neutron with one 
heavy-ion nucleus and also for optical potential model [21,22]. Although the non-relativistic 
Schrodinger equation with this potential has been solved for s-states [21] and the single- 
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particle motion in atomic nuclei has been explained quite well, the relativistic effects for a 
particle under the action of this potential are more important, especially for a strong-coupling 
system. Berkdemir et al [22] obtained the bound-state solution of the ID Schrodinger 
equation for the generalized WS potential by means of the NU method. We investigated the 
bound-state solutions of the ID KG equation with real and complex forms of the generalized 
WS potential [12]. Kennedy [23] obtained the scattering and bound-state solutions of the 
one-dimensional Dirac equation for the WS potential. Guo and Sheng [24] solved exactly the 
s(s)-wave Dirac equation {I — 0, i.e., K — —1 for spin and K — lior pseudospin symmetry) for 
a single particle with spin and pseudospin symmetry moving in a central WS potential. They 
obtained the energy spectra of the bound-states and the corresponding wave functions for the 
two-component spinor in terms of the hypergeometric functions. Alhaidari has developed 
a new two-component approach in the case of three-dimensional Dirac equation for the 
spherically symmetric potential and solved a class of shape-invariant Morse, Rosen-Morse, 
Eckart, Poschl-Teller and Scarf, potentials and given their rlativistic bound-state spectra 
and spinor wave functions [25-27]. Guo et al followed Alhaidari's approach and discussed 
the Dirac equation with WS and Hulthen potentials for spherical system and given their 
bound-state spectra and the spinor wave functions for s-states [28,29]. 

For the more realistic nuclear system where the nucleons are described in the relativistic 
mean field with the attractive S(/r) and repulsive vector Vi^r) potential, although some 
numerical techniques have been developed to Dirac equation, any analytic solution has not 
still been obtained for the Dirac- WS problem. In the special cases of spin symmetry (A = 
V — S = A =constant) and pseudospin symmetry {T^ = V + S = A =constant), Ginocchio 
et al solved triaxial, axial and spherical harmonic oscillators for the case A = and applied 
it to the study of antinucleons embedded in nuclei [30]. Lisboa et al studied the generalized 
relativistic harmonic oscillator for spin 1/2 particles and obtained the analytic solutions for 
bound states of the corresponding Dirac equations by setting A — [31]. Very recently, 
Ikhdair [32] studied the exact solution of the spatially dependent Dirac equation with the 
Rosen-Morse potential for arbitrary spin-orbit quantum number k. Under the conditions 
of the spin symmetry S ^ V and pseudospin symmetry S ~ — the bound state energy 
eigenvalues and corresponding upper- and lower-spinor wave functions are investigated in 
the framework of the NU method. Furthermore, We have solved the constant mass KG 
equation for the Eckart potential [32] and the spatially-dependent mass KG equation for the 
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Coulomb- like potential [33] and obtained the bound-state solutions of the energy eigenvalues 
and wavef unctions. 

In this work, we will solve the (3+1) dimensional Dirac equation for a particle trapped 
in the spherically symmetric generalized WS potential under the conditions of exact spin 
and pseudospin symmetry combined with approximation for the centrifugal and pseudo- 
centrifugal terms, and give the two-component spinor wavefunctions and the energy spectra 
for any arbitrary spin-orbit k bound states. The NU method [34] is used in the calculations. 
For the ID case, upon changing the values of the potential parameters from real to pure 
imaginary, we obtain Hamiltonians that may or may not be PT-symmetric. In addition, the 
spinor wavefunctions and the energy spectra of s-wave (/t = ±1) bound states and the non- 
relativistic limit are also discussed for the generalized and for a special case of the standard 
WS potential. 

The paper is structured as follows: In sect. 2, we outline the NU method and derive a 
parametric generalization version. Section 3 is devoted for the exact analytic bound state 
energy eigenvalues and two lower- and upper-spinor components wave functions of the Dirac 
equation with generalized WS quantum system obtained by means of the NU method. The 
spin symmetry and pseudo-spin symmetry solutions arc investigated using the NU method. 
In sect. 4, we study the cases k — ±1 {I — I — 0, s-wave) and the non-relativistic limit and 
compare with other wave equations and models. Finally, the relevant conclusions are given 
in sect. 5. 

II. THE NIKIFOROV-UVAROV METHOD 

The NU method has been used to solve the Schrodinger [19], KG [20,33] and Dirac 
[28,32] wave equations for central and non-central potentials. Let us briefly outline the basic 
concepts of the method [34]. This method was proposed to solve the second-order linear 
differential equation of the hypergeometric-type: 

a'^{z)u"{z) + a{z)T{z)u'{z) + a{z)u{z) - 0, (1) 

where the prime denotes the differentiation with respect to z, a{z) and 5(2;) are analytic 
polynomials, at most of second-degree, and r(s) is of a first-degree polynomial. Let us discuss 
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the exact particular solution of Eq. (1) by choosing 

U{z) = yn{z)(f)n{z), 



(2) 



resulting in a hypergeometric type equation of the form: 

a{z)y';{z) + r{z)y'^{z) + Xyn{z) = 0. (3) 

The first part yn{z) is the hypergeometric-type function whose polynomial solutions are 
given by the Rodrigues relation 

A (7" 

where is a normalization factor and p{z) is the weight function satisfying the condition 

[a{z)p{z)r^T{z)p{z), (5) 

with 

r{z) = t{z) + 27r{z),T'{z) < 0. (6) 

Since p{z) > and a{z) > 0, the derivative of t{z) has to be negative for bound states 
[32-34] which is the main essential condition for any choice of particular solution. The other 
part of the wave function is defined as a logarithmic derivative 

(l)'(z) 7t(z) 



4>{z) o-(^)' 
where 



<z) = \ {a\z) - t{z)\ ± {a\z) - r^z^ - a{z) + ka{z), (8) 
with 

k^X-7r'{z). (9) 

The determination of k is the key point in the calculation of 7t{z), for which the discriminant 
of the square root in the last equation is set to zero. This results in the polynomial Tr{z) 
which is dependent on the transformation function z{r). Also, the parameter A defined in 
Eq. (9) takes the form 

Xn^ -nT'{z) -]-n{n-l)a"{z), n = 0,l,2,---. (10) 
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At the end, the energy equation and consequently it's eigenvalues can be obtained by com- 
paring Eqs. (9) and (10). 

Let us now construct a parametric generalization of the NU method valid for any central 
and non-central exponential-type potential. Comparing the following generalized hypergeo- 
metric equation 

[z (1 - c^z)f u"{z) + [z{l- c^z) (ci - C2z)] u'{z) + (-Siz' + B2Z - B3) u(z) = 0, (11) 

with Eq. (1), we obtain 

t{z) = Ci - C2Z, a{z) = 2; (1 - csz) , a{z) = —Biz'^ + B2Z — S3, (12) 

where the parameters Cj and Bj (j = 1,2,3) are to be determined during the solution 
procedure. Thus, by following the method, we may obtain all the analytic polynomials and 
their relevant constants necessary for the solution of a radial wave equation. These analytic 
expressions are cited in Appendix A. 



III. SOLUTIONS OF THE DIRAC-GENERALIZED WS PROBLEM 

The Dirac equation for fermionic massive spin- 1/2 particles moving in an attractive scalar 
potential S{r) and a repulsive vector potential V{r) can be written as [35] 

[ca • p + /3 (mc" + S{r)) + V{r) - E] ^{r) = 0, V(r) = V(r, ^, 0), (13) 

where E is the relativistic energy of the system, p = —ihV is the momentum operator, and 
a and (3 are 4x4 Dirac matrices 

« = I ' h = I ^ ° I , ^^1= I ° M , (72= ° ~ I , (73- I ^ ° I . (14) 





where I denotes the 2x2 identity matrix and cr are three-vector Pauli spin matrices. For 
a spherical symmetrical nuclei, the total angular momentum operator of the nuclei J and 
spin-orbit matrix operator K = — /3 (cr • L -|- I) commute with the Dirac Hamiltonian, where 
L is the orbital angular momentum operator. The spinor wavefunctions can be classified 
according to the radial quantum number n and the spin-orbit quantum number k and can 
be written using the Pauli-Dirac representation: 

I friK \ 1 I ^nK{f)Yjj^i^{r) 



where Fni^{r) and GnK{f) are the radial wave functions of the upper- and lower-spinor 
components, respectively, Y^^^^ij) and Yjm(-K)^) denote the spin spherical harmonic func- 
tions coupled to the total angular momentum j and it's projection m on the z axis and 
/(/ + 1) = 1) and /(/ + 1) = K (k — 1). The quantum number n is related to the 

quantum numbers for spin symmetry / and pseudospin symmetry / as 



K — 



(/ + !) = -(j + l), (si/2,P3/2, etc.), J=/ + |, aligned spin (k < 0) , 
+/ = + (i + \) , (pi/2, c?3/2, etc.), j = l-\, unaligned spin {k > 0) 



(16) 



and the quasi-degenerate doublet structure can be expressed in terms of a pseudospin angular 
momentum 's — 1/2 and pseudo-orbital angular momentum I which is defined as 



K 



-^ = -(j + l), {si/2,P3/2, etc.), j ^1-1/2, aligned spin < 0) , 

/ + l) = + (j + i) , (Pi/2, d3/2, etc.), j = T+ 1/2, unaligned spin {k > 0) . 

(17) 

where k — ±1,±2,--- . For example, (3si/2! 2^3/2) and (3pi/2, 2^3/2) can be considered as 
pseudospin doublets. 

Substituting Eq. (15) into Eq. (13), we obtain two radial coupled Dirac equations for 
the spinor components 



+ 7) ^nnir) = [mc^ + E„« - A(r)] G„«(r), 



d K 



Gn.{r) = [mc' - En. + E(r)] F„«(r), 



(18a) 



(18b) 



where A(r) = V{r) — S{r) and S(r) = V{r) + S{r) are the difference and sum potentials, 
respectively. By eliminating Gnnir) in Eq. (18a) and Fn^if^) in Eq- (18b), one is able to 
obtain two second-order differential equations for the upper- and lower-spinor components 
as follows: 



d^ k{k + 1) 1 



d"^ _k(k-1) 1 
dr'^ ^2 



C/_(r)C/+(r) 



dMr) (d_ I k\ 
dr \dr r ) 

mc2 + E^^ - A(r) 



Fn.{r) = 0, F„«(0) = 0, 
(19a) 

Gnn{r) = 0, ^^.(O) = 0, 

(19b) 

where U-{r) — mc^ -\- En. — A(r) and C/+(r) = mc^ — En. -|- E(r), are the difference and 
the sum functions, respectively. Prom the above equations, the energy eigenvalues depend 
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h'6 



U4r)U+{r)+ '^^ ^''^ '■^ 



mc2 - Enn + S(r) 



on the quantum numbers n and and also the pseudo-orbital angular quantum number I 
according to k{k — 1) = l{l + 1), which imphes that j — I ± 1/2 are degenerate for I ^ 0. 
The above non-linear radial wave equations having very complicated solutions are required 
to satisfy the necessary boundary conditions F„k(0) = ^^^(O) and F„«;(oo) = G„«;(oo) for 
bound state solutions. 

In this context, we take the sum potential in the form of an attractive generalized WS 
potential, i.e., E(r) = yawsif) [36]. The interaction among nuclei is commonly described by 
using a potential which consists of the Coulomb and the nuclear potentials. It is usually taken 
in the form of WS potential. Here we take the following Hermitian real- valued generalized 
WS potential which is specified by the shape (deformation) parameter, g, [13,22,37] 

Vcwsix) = Vg{x) = -Vo ^ Ro/a = a, (r - Rq) /Rq = x, i?o > a, g > (20) 

where r e (0, oo) or x e (—1, oo) refers to the center-of-mass distance between the projectile 
and the target nuclei. The relevant parameters of the inter-nuclear potential are given as 

1/3 

follows: Rq = tqAq is to define the confinement barrier position value of the corresponding 
spherical nucleus or the range of the potential well, Aq is the atomic mass number of target 
nucleus, Tq is the radius parameter, the parameter Vq is the potential depth, a is the surface 
thickness and has to control it's slope, which is usually adjusted to the experimental values 
of ionization energies. Note further, q is a real shape (deformation) parameter, the strength 
of the exponential part other than unity, set to determine the shape of potential and is 
arbitrarily taken to be a real constant within the potential. In addition, it should be noted 
that the spatial coordinates in the potential are not deformed and thus the potential still 
remains spherical. 

It is worth noting that under radial coordinate transformation, r — > r -|- A, then the 
generahzed WS potential in Eq. (20) changes into the standard WS potential (when q is 
taken equal to 1 in the calculation) but with the displacement parameter A satisfies the 
expression cxp(A/a) = q and with a field strength Vq = VoCxp(— A/a) [13,22]. The sense 
of generalization or deformation of the potential becomes clear. For completeness, it could 
be stated that if A is positive (corresponding to g > 1) then one may need to impose the 
condition on the choice of A, that is, |A| <^ Rq. 

Obviously, for some specific q values this potential reduces to the well-known types, such 
as for g = to the exponential potential and for g = — 1 and a — to the generalized 
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Hulthen potential (cf. [11,13,22] and the references therein). Obviously, the solutions in 
Ref. [22] are at best vahd for Rq — 0, in which the potential can be expanded in terms of 
hyperbohc functions [19,22]. The standard WS potential turns to become the well-known 

Rosen-Morse potential shifted by the term — Vo/2 (cf. Ref. [22] and the references therein), 
that is, VsRAiir) = —Visec h'^{ax) + V2tanh(Q;a;) — V3, where Vi = C/4 and V2 = V3 = Vq/2 
[32]. 



A. Spin symmetric solution 



In the case of exact spin symmetry S{r) ~ V{r) {dA{r)/dr — 0, i.e., A(r) 
constant), Eq. (19a) can be approximately written as 

where + — / (/ -|- 1) , n = / for «; < and k = — (/ -|- 1) for k > 0. The spin 
symmetric energy eigenvalues depend on n and k, i.e., E^k — E{n, k {k + 1)). For I ^ 0, 



^ [miry + K« - A] [miry - E^,, + E(r)] | F„«(r) = 0, (21) 



the states with j — I ±1/2 are degenerate. This is the exact spin symmetry. Taking the 
E(r) = 2y(r) — > Vowsif) ^ mentioned in Ref. [36] enables one to reduce the resulting 
relativistic solutions into their non-relativistic limit under appropriate transformations. We 
are set out to obtain bound state solutions (relativistic energy spectrum and upper- and lower 
spinor wavcfunctions) of a spin- zero particle for a four parameter {Vo.q, a, Rq} generalized 
WS potential by means of the NU method. Moreover, if k is not too large, the case of 
the vibrations of small amplitude about the minimum, we can then use the approximate 
expansion of the centrifugal potential near the minimum point r — Rq ix — 0) as [38] 

Kin + l) K (k -f- 1) 



i?0 



where 



Rq 



— exp(— axj 
1 -|- qexpi—ax) 



[1 + xy 



+ D2 



— exp(— 



1 -|- qexpi—ax) 



Di = 2 [exp(Q;i?o) + 1] |3 
D2 = [exp(Q;i?o) + 1]^ 



exp(— Q;i?o) 
aRo 

1 + exp(— Q;i?o) 
aRo 

1 -|- exp(— Q;i?o) 
aRn 



iaRo 



1 + exp{—aRo) 
- (3 + aRo) 



3 -|- aRo — 



3 — aR 



1 + exp(— Q;i?o) 
aRo 

2aRo 



1 -|- exp(— Q;i?o) 



(22) 

(23a) 
(23b) 
(23c) 



and higher order terms are neglected. It is worth noting that for k 7^ case, we have to use 
an approximation for the centrifugal term similar to the non-relativistic cases which is valid 
only ioT q — 1 value [6,38]. However, for s- waves, we remark that the problem can be solved 
exactly and the solution is valid for any deformation parameter q. 

We define the following new dimensionless parameter, z{x) = — e^"^ G [— e^°/",0], which 
maintains the finiteness of the transformed wave functions on the boundary conditions. 
Thus, using Eqs. (22) and (23) , we can reduce Eq. (21) to generalized equation of the 
hypergeometric type for the upper-spinor component F„«;(r), 



(1 - qz) d 1 t n 1 n 2 \ 



where -^^^(r) = F^^XA introduce the definitions 



F^,{z) = 0, F„«(0) = F„«(oo) = 0, 

(24) 



-UK 



\m^c^ -El^-A {mc^ - E„,)] + -^D^ > 0, (25a) 



Pi - <fel^ - ^ (^^' + E^^-A)-^ {qD, - D,) , (25b) 
/52 = 2qel^ - ^ (mc^ + K« -A)- (25c) 

il C Tq 

with uj = k{k + 1) where k = ±1, ±2, ■ ■ ■ , for bound states (i.e., real enn)- Before we can 
proceed, it is necessary to compare the last equation with Eq. (1) to obtain the following 
polynomials: 

T{z) = l-qz, a{z)=z{l-qz), a{z) = -(5^z^ + (32Z - el^. (26) 

We follow Appendix A to calculate the specific values of the parametric constants and then 
display them in Table 1 for the present potential model. Also, with the aid of Table 1, the 
key polynomials given in Appendix A take the following particular analytic forms: 

t:{z) = - I (1 + 2e^^ + ^ z, (27) 

k = p2-q i^Enn + £nK, (28) 

t{z) = 1 + 2£„, - ? (2 + 2£„, + z, (29) 
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where t'{z) = -g (2 + 2£„« + ^ < with ^ = ^1 -r 



^'^"'i^2. We insert the values of the 

i 

constants given in Table 1 into the energy equation cited in Appendix A and then obtain 

2 



+ 



2a 



1 + 2n + W 1 + ^^D2 



(30) 



where k = ±1, ±2, • • • . Hence, the above equation gives explicitly the energy equation with 
exact spin symmetry for arbitrary spin-orbit coupling quantum k of the Dirac equation as 
follows 



Rl 



-Dn 



ahc 



-1 2 



l + 2n+^l + ^D,) 



(31) 

The energy level Enn is determined by energy equation (31), which is a rather complicated 
transcendental equation. Now, let us consider a few special cases of much concern, (i) If we 
choose q — 1, the potential (20) turns to the shifted WS potential: 



V{x) = -V^o + Y 



+ e- 



(32) 



and then it's energy spectra yield 



I t2 2 2 

+ Ti c a 



{mc^ - Er,.) [mc^ + E^,. - A) 



^{mc^ + En,-A) + ^(D,-D2) 1 



-Dn 



(l + 2n+^l + ^D, 



4cja2 
l + ^r7^/^2 



n 2 



Rq 



(33) 



(ii) If we choose q — —l, the potential (20) turns to the standard shifted Hulthen potential: 

Vo 



and then the resulting energy eigenvalues become 

2 

(mc^ - Enn) {mc^ + E^n - A) ^ 52^ ^0 



(34) 



I i;2 2 2 

+ rica 



^ (mc^ + K« - A) + 1^ {Di + D2) 1 



-1 2 



(l + 2n+^l + ^Z^2 



^^,^l + 2n+yi + ^D. 



(35) 
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(iii) If we choose g — > 0, the potential (20) turns to the exponential potential: 

V{x) = -Foe-"", (36) 

the eigenvalues expression (31) does not give an explicit form, i.e., the NU method is not 
applicable to the exponential potential (36). Note that for this potential there is no explicit 
form of the energy expression of bound states for Schrodinger [9] , KG [12] and also Dirac 
[8] equations. 

In addition, for the s-wave {k — —1) and V{r) — S{r) (i.e., A — 0), we obtain 

2 

, n = 0,l,2,--- , (37) 



*2 2 2 
2 4 ■^^2 lb ^ a 

m c — Et — 



\/o (mc2 + ^ g(n + r 



2ri^c^ n + 1 2a2 
and it can be seen easily that while the field strength Vq 0, the energy states yield: 

= d=^^4a2m2c4 - n^c^ {n + 1)', n = 0,l,2,---, (38) 

for particles and anti-particles. Note that in the above equation there exist bound states for 
the ground and excited states in — 0, 1) which are — ■^\/?>m(? /2 and Ei — 0, respectively, 
for positive q values and a = Ac, where Ac = h/mc denotes the Compton wavelength of the 
Dirac particle. Otherwise, there are no bound states for n > 2 states. 

On the other hand, for the same value of a and negative q values when Vq ^ 0, all energy 
eigenvalues go to zero. If the value of q is increasing, all positive bound states go to zero, 
from (38), asymptotically. 

An inspection of the energy expression given by Eq. (37), for any given a, shows that we 
deal with a family of generalized WS potentials. The sign of Vq does not effect the bound 
states. The spectrum consists of complex eigenvalues depending on q. As we shall see the 
role played by the range parameter a is very crucial in this regard. Of course, it is clear 
that by imposing appropriate changes in the parameters {a,Vo,q} , the energy spectrum in 
Eq. (37) for any modified parameter can be also calculated by resolving Dirac equation for 
every parameter change. 

The uppcr-spinor wave functions for Fnnir) will be presented. In order to establish FnK.{r), 
use will be made of Appendix A and Table 1. Hence, the first part of wave functions reads: 

<P^{z) = z'-{l - g^)^(i+«), Sn. > 0, e > -1- (39) 
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In addition, to find the function, yniz), which is the polynomial solution of hypergeometric- 
type equation, we firstly calculate the weight function: 



and thus the second part of wave functions (7) can be obtained as 



(40) 



yn{z)=DnZ-^'-'^{l-qz)-^^^ 



1-z 



, e>-l, 



(41) 



where Dn is a normalization constant. In the limit q ^ 1, the polynomial solutions of 
yn{z) arc expressed in terms of Jacobi Polynomials, which is one of the classical orthogonal 
polynomials, with weight function given by Eq. (40) for z G [0, 1], giving yn{z) ~ P^^"'^'^\l — 
2z), 2eni^, ^ > — 1 [39]. Thus the associated uppercomponent Fni^{z) for arbitrary the spin- 
orbit coupling quantum number k can be obtained by substituting Eqs. (39) and (41) into 
Eq. (2) as 

Fn.{z) = Mnnz'-il - g^)5(i+«)pi2£„.,0(i _ 2qz) 
= AA„„^""-(l-g^)^ (1+^)2^1 (-n,l + 2£„„ + e + n;l + 2£„„;g^), (42) 

where z{r) = — e"^''"'^)/" and jV„k are normalization constants calculated in Appendix B. 

Before presenting the corresponding lower-component GnK{i^), let us recall a recurrence 
relation of hypergeometric function, which is used to solve Eq. (18a) and present the corre- 
sponding lower component GnK{i^), 



d 
dz 



2F1 {a;b;c;z) 



— ) 2F1 (a + l;6+l;c+l;^). 



(43) 



with which the corresponding lower component GnK{f^) can be given by solving Eq. (18a) as 
follows 



mc2 - A 



dr r 



+ - 



2a (1 + ge-('-'f^)/«) r 



X 



2F1 (-n, 1 + 2enK + ^ + n; 2£„«, + 1; -qe 



-{r-Ro)/a 



) 



qn 



;i + 2£„« + ^ + n) (_e-(^-«°)/«)'"'^'^' (1 + ge-('-«°)/«)^^'+^^ 



a {2enK + 1) {mc^ + E^n - A) 
X (-n + 1; n + e + 2 (1 + £„«) ; 2 (1 + £„,) ; -qe'^^-^)!'') . (44) 
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Here, it should be noted that the hypergeometric series 
2F1 {—n, 1 + 2enK + ^ + 2enK + 1; —qe'^^~^^^"') does not terminate for n = and 
thus does not diverge for all values of real parameters ^ and 

For A > mc^ + Enn and < mc^ or A < mc^ + E^^ and E^^ > mc^, we note that 
parameters given in Eq. (25a) turn to be imaginary, i.e., < in the s-statc {k, = — 1). 
As a result, the condition of existing bound states are > and ^ > 0, that is to say, in 
the case of A > mc^ + EnK and < mc^, bound-states do not exist for some quantum 
number k such as the s-state {k — —1). Of course, if these conditions are satisfied for existing 
bound-states, the energy equation and wave functions are the same as these given in Eq. 
(31) and Eqs. (42)-(44). 



B. Pseudospin symmetric solution 

Under the condition of the pseudospin symmetry S{r) ~ —V{r) (i.e., dTi{r)/dr — 0, or 
E(r) — A' — constant), Eq. (19b) can be exactly written as 

{ ^ - ^ - ^ [^c' + - ^(0] [^c' - + ^] } Gn.{r) = 0, (45) 

where u) — k{k — 1) — l{l -\- 1), the energy eigenvalues E^^ depend only on n and I, i.e., 
EnK — E{n,l{l + 1)). Taking the A(r) = Vawsif) allows us to reduce our results to the 
non-relativistic limit. For I ^ 0, the states with j — I ±1/2 are degenerate. This is the 
exact pseudospin symmetry. We follow the procedures in the previous subsection to obtain 
Dirac equation satisfying GnK(^); 

1 — qz) d 1 



^ ^ W^)d~z + [.(1 - qz)f ("^^"' + " ^-^), 



C(z) = 0, C(0) = G'„,(oo) = 0, 

(46) 

where GnK{f) — Gnniz) and we have used the definitions 



^nn = \ hTl [^'c' - + A' {mc^ + En^)] + -^Do > 0, « = ±1, ±2, • • • , (47a) 
A = I'eL + ^ {rnc' - E^. + A')-^ {qD, - D,) , (47b) 

lb C -Kg 

(32 = 2qel, + {mc' - E^. + A') - —^D,. (47c) 

ll C -ftQ 



14 



To avoid repetition in the solution of Eq. (46), a first inspection for the relationship be- 
tween the present set of parameters {suk, Pi, P2) and the previous set {snK, Pi, P2) tells us 
that the negative energy solution for pseudospin symmetry, where S{r) ~ —V{r), can be 
obtained directly from those of the positive energy solution above for spin symmetry using 
the following parameter mapping [39-41]: 

Fn.{r) ^ Gn.{r), V{r) ^ -V{r) (or Vo ^ -H), ^n. ^ -^n. and A ^ -A'. (48) 

Following the previous results with the above transformations, we finally arrive at the energy 
equation. The relativistic transcendental energy equation is 

[m'c' - El + A' {mc' + E^^)] = -^Do 



Rq 



( ahc^ ^ 



+ - 



V 1 



.^^[m^-E^. + A') + ^{qD,-D,) g / ^ n 



(1 + 2n + ^1 + f§D^) \^ y q^l ' 

(49) 

and the lower-spinor wave functions 

= K^z^-'^il - qz)'^('^^ ^F, (-n, 1 + 2?'„^ + ^ + n; 2£„« + 1; qz^ , (50) 



with ^ = -^1 + ^^'d2 and is defined in Eq. (47a). 

IV. DISCUSSIONS 

Now, let us study three special cases. We first study the s(s)-states {1 = 1 = 0, i.e., 
K = =1=1 ). In this case, wc have the spin-orbit coupling term k,(k + l)/r^ = 0, and also the 
corresponding approximation to it in Eq. (22). The corresponding energy equation reduces 
to the s-states {k = —1), i.e., 

2 



aVo {mc^ + - A) 1. . 

qh'c^l + n) ^a^'^''^ 



(51) 

and the upper-spinor component of the wave functions: 

Fn{z) = Kz'-{1 - qz)Pi^'-^'\l - 2qz) 
15 



Nnz'-'il - qz) 2F1 (-n, 2(1 + £„) + n; 1 + 2£„; 5^) , 



(52) 



with 



£n = ^^m^c^ - - ^ (mc2 - K(-i)) > 0, 



(53) 



where jV„ is calculated in Appendix B. As mentioned above, in the s-wave {k = —1) the 

condition of existing bound-states is for A < mc^ + En(-i) and -E'„(-i) < mc^. Furthermore, in 
the nonrelativistic limit with the mapping mc^ — En{-\) —E^q and (mc^ + £"„(_!)) 
1^, then we have 



2^2 



En — 



h a 
""8/7 



2/i H ^ (1 + n) 



hUil + n) ' a2 J ' ^^^^ 
Second, we study the special case A(r) = ^4 = 0. If so we have V{r) — S{r) — l/2E(r) 
and it turns to the KG solution. Obviously, in this case the energy equation given in Eq. 
(31) reduces to the energy equation of arbitrary k state Dirac equation for equal scalar and 
vector WS potential as follows 



2 4 

m c 



Rl 



( ahc^ ^ 



+ - 



V 1 



^ w + £;„,) + ^ (gL*! - ^2) 



-1 2 



he 



l + 2n + Jl + ^D, 



, (66) 



and the upper component of the wave functions 



XnnZ'"'{l - g^)^(^+^°) 2F1 {-n,l + 2£n« + + 2£n« + l', Qz) , 



^ {m^c^ - El,) + 



rir, 



1^ , 4^(/^ + l)a^ 



where E^, < m^c^ + '^^g^ Dq is the essential condition for existing bound-states. 
Third, the non-relativistic energy state limit for arbitrary / state are 

_ h'lji + 1) 



(56) 
(57) 



2^2 



h a 



I ^ Alil + l)a?^ 



(58) 
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The above result is identical to Eq. (23) in Ref. [22] where Berkdemir et al. used the usual 
approximation to the centrifugal term in the potential expression (10) (cf. [42] and Eq. (2) 
in J. Math. Chem. 42, 461 (2007)). 

It is worthy to note that in the calculations of Ref. [22], Rq was neglected. Accordingly, 
the solutions of the energy spectra Eq. (23) of the original paper [22] are at best valid for 
i?o = in which case the standard WS potential Eq. (10) in the original paper reduces to 
the shifted Rosen-Morse (RM) potential (cf. Eq. (3) in Phys. Rev. C 74, 039902(E) (2006)). 
The additional potential besides the standard Ws potential considered by Berkdemir et al. 
[22] provides the flexibility to construct the surface structure of the related nucleus [42]. 
Thus, the non-relativistic solutions obtained in [22] are only reasonable for the hyperbolic 
[43] exponential (RM) potential [32] , not WS potential. This clear when we rewrite Eq. (20) 
in the following form 

VGws{r) = -V' ' a = -, = ^oe^"/", q = e^^'\ r e (0, oo), (59) 

when i?o = 0, it implies that g = 1 and then the above potential reduces to the standard WS- 
type potential. In addition, the authors of Ref. [22] approximated the centrifugal potential 
term ^ l{l + iW ^^^'IZrf [13,43], where C ^ l{l + l)a^ [20]. However, in the present 
work, Eq. (58) contains the width of the potential Rq. Also, an expansion for the centrifugal 
potential term has been performed around the point r ^ Rq {x — 0) [38] , and without loss 
of generality we put x = {r — Rq) /Rq at the end of our calculations. 

The empirical values found by Percy et al. are given as Rq = 1.285 fm and a = 0.65 
fm [44]. In addition, the following WS potential strength parameter is Vq ~ 40.5 + 0.13/lo 
MeV in the non-relativistic limit. Here, Ao is the atomic mass number of target nucleus and 
is defined through Rq — tqAI/^. On the other hand, the associated upper-spinor component 
of the wave functions is 

= Nr,z'-^ (1 - g^) (-n, 1 + 2eni + Ci + n; 2eni + l;qz), (60) 

where 

m + i) I m + i)a2 ^ ^ 

e„i = aJ-^Er,i + ^^Do>0, Ci^Jl+ ^^.^^ D^, Z = 0,l,2,---, (61) 
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where Eni < 2hr^ -^o essential condition for existing bound-states. As a numerical 

example, we impose appropriate values for the parameters in Eq. (51) to calculate the bound 
state energies of the spin symmetry generalized WS potential for special case k, = —1 and 
using h = c = 1. The results obtained by using the following parameters Vq = 2.2, m = 15, 
A — —5 and a — 1.425 are given in Table 2. The condition of existing bound states is the 
A < 0. When the A > 0, there are no bound states in the hmit of exact spin symmetry. For 
g = ±1, there is only one attractive bound state Eo_i — (A + m) . 



V. CONCLUSIONS 



We have discussed the approximate bound state solutions of the Dirac equation for the 
generalized WS potential with any arbitrary spin-orbit k state under the conditions of the 
spin (pseudospin) symmetry V — S — A {V + S — A) hy means of the NU method com- 
bined with the approximation for the centrifugal term. By setting V + S {V — S) to the 
spherically symmetric WS potential, we have derived the solutions of the Dirac equation for 
the relativistic energy eigenvalues and associated two-component spinor wave functions for 
arbitrary spin-orbit n state that provides an approximate solution to the spin- and pseudo- 
spin symmetric Dirac equations. The most stringent interesting result is that the present 
spin and pseudo-spin symmetric can be easily reduced to the KG solution once S{r) — V{r) 
and S{r) — —V{r) {i.e., A — 0), respectively. The non-relativistic hmits of our solution 
are obtained by imposing appropriate transformations. The resulting solutions of the wave 
functions are being expressed in terms of the Jacobi polynomials. If we choose the spin-orbit 
quantum number k = —1 {k, = 1) for spin (pseudospin) symmetry, the problem reduces to 
the exact s(?)-wave Dirac solution. The generalized Hulthen potential bound state solutions 
are simply derived when letting q ^ —q. We have also discussed the relation between the 
non-relativistic and relativistic solutions and the possibility of existing the bound states. 
It should be noted that the numerical calculation for energy levels involved in Eq. (51) 
is terribly sensitive to the choice of those parameters. In Table 2, we choose _^ as the 
physical solution for the transcendental equation (51). 

At the end, the solutions that constitute the main results regarding the energy equations 
(31) and (49) for the spin and pseudospin symmetry, respectively, may have some interesting 
applications in many areas in physics. For example, the work is helpful to understand 



18 



spectroscopy with high field physics and useful to understand the nuclear properties like 
nuclear scattering systems [45]. In addition, the present results play an essential role in 
microscopic physics, since it can be used to describe the interaction of a nucleon with a 
heavy nucleus [21,22]. In the non-relativistic limits, the energy eigenvalues, Eq. (58), is 
physical and is in good agreement with the results obtained previously by other methods 
and works [22]. 

Acknowledgments 

We wish to thank the anonymous referees for their invaluable suggestions and comments 
that have improved the paper greatly. We are also grateful for the partial support provided 
by the Scientific and Technological Research Council of Turkey (TUBITAK). 

APPENDIX A: PARAMETRIC GENERALIZATION VERSION OF THE NU 
METHOD 

We complement the theoretical formulation of the NU method in presenting the essential 
polynomials, energy equation and wave functions together with their relevant constants as 
follows [46]. 

(i) The key polynomials: 

Tt{z) = C4 + C5Z - [{^/c^ + C3^/C^) z - ^/c^ , (Al) 

= - (c7 + 2C3C8) - 2V^. (A2) 

t{z) = 1 - (c2 - 2C5) z-2 [(v^ + C3V^) z-^, (A3) 

t\z) = -2c3 - 2 (v^ + cs V^) < 0, (A4) 

(ii) The energy equation: 

(C2 - C3) n + 03^2 - (2n + 1) C5 + (2n + 1) + €3^) + C7 + 2C3C8 + 2^/^ = 0. (A5) 

(iii) The wave functions: 

p{z)=z''''^{l-cszy'\ (A6) 

cf>{z) = z'^^'il - cszY'^ C12 > 0, ci3 > 0, (A7) 
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Vr^iz) = P^-io.'^")(i _ 2csz), Cio > -1, Cii > -1, (A8) 

U{z) = MnZ^^^{l - C3Zy''Pj^'''°'^''\l - 2C3Z), (A9) 

where Pn^''^\x), ji > > —1 and x e [—1,1] are the Jacobi polynomials with 

P^",/3)(l _ 2s) = i^lik 1 + q; + /3 + n; a + 1; s) , (AlO) 

77/. 

and is a normalization constants. Also, the above wave functions can be expressed in 
terms of the hypergeometric function as 

Unniz) = NnnZ^^^il - C^zf^^ ^Fi {-n, 1 + Cio + Cn + n; Cio + 1; C3Z) , (All) 

where Cu > 0, C13 > and z G [0, 1/C3] . 
(iv) The relevant constants: 

C4 = ^ (1 - Ci) , C5 = ^ (C2 - 2C3) , C6 = C5 + Bi, 

C7 = 2C4C5 - B2, Cg = C4 + ^3, Cg = C3 (C7 + C3C8) + Ce, 

2 

Clo = ci + 2c4 + 2a/c^ - 1 > -1, Cn = 1 - ci - 2c4 H a/c^ > -1, 

ca 

C12 = C4 + > 0, ci3 = -C4 + — (v^ - C5) > 0. (A12) 



APPENDIX B: NORMALIZATION OF THE RADIAL WAVE FUNCTION 

In order to find the normalization constants Muk, we start by writting the normalization 
condition: 

(^L [ z^'---\l - z)^+^ [P^'^"-«(l - 2z)]^dz = 1. (Bl) 
Jo 

where q = 1- Unfortunately, there is no formula available to calculate this key integration. 
Neveretheless, we can find the explicit normalization constant A/"„k- For this purpose, it is 
not difficult to obtain the results of the above integral by using the following formulas [46,47] 

J (1 - s s"-^ 2F1 {a, P; 7; as) dz = ^^^^^^^ 3F2 {ly, a, + z/; 7; a) , > -1, z/ > -1, 

^ (B2) 
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and 2F1 (a, b; c; z) = E^o ^^^W^f- Following Ref. [46], we calculate the nor- 

malization constants: 



where 



ar(2e„, + l)r(e + 2) ^ (-1)^ (1 + n + 2£„. + Q)^ r(n + m) ^ 



r(n) 



^ m! (m + 2£„«) !r (m + 2£„« + C + 2) 



-1/2 



(B3) 



Cn,^ = (2£n« + m, -n, n + 1 + C + 2£„„; m + 2£„« + ^ + 2; 1 + 2£„«; 1) , (B4) 



Furthermore, the normalization constants for the s-wave can be also found as 



mathcalNnK — 



where 



ar{2enn + l)r(e + 2) ^^ (-1)" (1 + n + 2e„. + 0)„ T{n + m) 

2^ z rrrr^ z : r^O„ 



r n 



m=0 



m! (m + 2£„«)!r (m + 2£„« + C + 2) 



-1/2 



(B5) 



C„ = 3F2 (2£n + m, -n, n + 2£„ + 2; m + 2£„ + 3; 1 + 2£„; 1) , (B6) 
and 6n is given in Eq. (53). 
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TABLE I: The specific values for the parametric constants necessary for the present potential. 



Constant Value Constant Value 



Ci 


1 


C2 




C3 


Q 


C4 





Ct; 
^0 


Q 
2 


«^ 


1 fo2 + 4/3i) 




-(32 


C8 




eg 




ClO 




Cll 




C12 




Cl3 






B2 


02 


^3 


£2 



TABLE IL The bound state energy levels E^k for the special case k = —1. 





q = l 






q = 2 






q = -l 






q = -2 






n 






1 






-1 


Efi-i 




-1 






1 





-10.197602 


-19, 


.996367 


1.337420 


-19, 


.996426 


-9.561001 


-19, 


.996589 


2.018956 


-19, 


.996536 


1 


0.985745 


-19, 


.985463 


9.849615 


-19, 


.985698 


2.349482 


-19, 


.986350 


10.803623 


-19, 


.986141 


2 


6.597152 


-19, 


.967274 


12.216763 


-19, 


.967804 


8.327343 


-19, 


.969273 


13.247481 


-19, 


.968802 


3 


9.328912 


-19, 


.941777 


13.119033 


-19 


.942722 


11.239603 


-19 


.945338 


14.180142 


-19, 


.944500 


4 


10.775321 


-19, 


.908939 


13.534261 


-19, 


.910418 


12.784029 


-19, 


.914518 


14.610668 


-19, 


,913204 


5 


11.602814 


-19, 


.868715 


13.742899 


-19, 


.870854 


13.670315 


-19, 


.876777 


14.828469 


-19, 


.874879 


6 


12.102341 


-19, 


.821054 


13.847458 


-19, 


.823977 


14.208309 


-19, 


.832071 


14.939340 


-19, 


.829478 


7 


12.412538 


-19, 


.765892 


13.892224 


-19, 


.769728 


14.545655 


-19, 


.780348 


14.988985 


-19, 


.776946 


8 


12.605105 


-19, 


.703154 


13.898750 


-19 


.708035 


14.758727 


-19 


.721546 


14.999654 


-19, 


.717219 


9 


12.719786 


-19, 


.632753 


13.878495 


-19 


.638818 


14.889869 


-19, 


.655595 


14.983185 


-19, 


.650223 


10 12.779911 


-19, 


.554593 


13.837978 


-19, 


.561982 


14.963949 


-19, 


.582415 


14.946314 


-19, 


.575874 
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